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, Abstract. This course consists of two lectures. In the first lecture I discuss why a non perturbative 

formulation of QCD is needed ,and I show that lattice formulation copes with this need, even if it 
mainly produces numerical results. In the second lecture I discuss how lattice can help to understand 
the deconfinement transition. Such understanding is also important to predict parameters that can 
5-H ■ help in the interpertation of heavy ions high energy experiments. 

Q-i' 

<: 

<r> [ INTRODUCTION 

>• [ Vacuum is by definition the ground state of a field system:it is stable against quantum 

fluctuations. 

An exact knowledge of the ground state provides complete information on the sys- 
tem.This property goes under the name of reconstruction theorem,the exact statement 
^ \ being that from the field correlators 

(0\T(<p 1 (x 1 )...( Pn (x n )\0) (1) 

a . 

the Hilbert space and the matrix elements of physical observables can be constructed. 
\ Textbook quantization is perturbative .The Lagrangean is split in a free term Jz? and 

^3 ■ a perturbation Jz? 7 

>: ' j? = j? +j?j (2) 

Jf is exactly solvable ,and defines as vacuum the ground state of Fock space,Le. empty 
space. ££ l is a small perturbation producing scattering between fundamental particles, 
and small changes in the wave function of the ground state. 

In QCD J?? describes free gluons and quarks and their interactions. For sure this is 
not a good starting point , since quarks and gluons are confined.Fock vacuum is not the 
ground state,and therefore it is not stable against perturbations. 

This is most probably the reason why the perturbative expansion is not conver- 
gent,not even as an asymptotic series,in spite of the fact that it seems to work at small 
distances [[I]]. 

The knowledge of the true vacuum is,in particular,relevant to the understanding of 
one of the most intriguing properties of QCD, namely confinement of color. 

In spite of the well established evidence that quarks and gluons are the fundamental 
constituents of hadronic matter, they have never been observed in Nature as free parti- 
cles. 

All particles produced in particle reactions are color singlets. This property is known 
as confinement of color. 



Since the pioneering papers of Gellmann in which quarks were introduced as funda- 
mental constituents of hadrons ,quarks have been searched for in particle reactions and 
in Nature.None has been found,so that experiments establish upper limits, which can be 
found in Particle Data Group[Q]. The cross section o q for the inclusive reaction 

p + p^q(q)+X (3) 

in which a q, (q) is detected as a fractionally charged particle has an upper limit 

a q < l(T 40 cm 2 (4) 

to be compared to the total cross section o T at the same energy o T ~ 10 25 cm 2 . 
In the absence of confinement o q /o T would be of the order of unity,and is instead 
o q /o T ~ 1(T 15 . 

Similarity the abundance n q of free quarks in Nature is measured by looking for 
particles of fractional charge q = ±1/3, ±2/3 in Millikan like experiments. 

No quarks have been found,and the resulting upper limit is n q /n p < 10~ 27 corre- 
sponding to the analysis of ~ Ig of matter. 

In the absence of confinement the Standard Cosmological Model (SCM) predicts 
n q /n p ~ 10 _12 []3p. Again a factor 10~ 15 between upper limit and expectation. 

Such a small factor cannot have a natural theoretical explanation ,except if the actual 
value of the above quantities is exactly zero. Confinement is an absolute property,and this 
can only be explained in terms of a symmetry property of the vacuum. The knowledge of 
the true vacuum is necessary to understand confinement .A quantization procedure that 
is not based on perturbation theory is needed. 



FEYNMAN PATH INTEGRAL [f|] AND LATTICE FORMULATION 

OF QCD[§. 

A quantum system is defined by the canonical variables ,q,p and the Hamiltonian H. 
Solving the system means to construct a Hilbert Space on which q,p act as operators 
obeying the equations of motion and the canonical equal-time commutation relations. 

[<liiPj} = i5 ij (5) 

A ground state must exist. 

Field theory is a special case. The fields <p ; (x, t) play the role of the q's ,their conjugate 
momenta Tl^xj) = dJf/d(d <p i (x,t)) the role of the p's. The system has infinitely 
many degrees of freedom. To simplify the notation we shall refer in what follows to a 
system with one q and one p , since our arguments will apply to any number of degrees 
of freedom. Problems with divergencies ,arising from infinite number of degrees of 
freedom are technical in nature and do not affect the arguments below.The same is 
for problems in defining the conjugate momenta, which are a consequence of gauge 
invariance. 

Since the \q)'s are a complete set of states ,the knowledge of the amplitudes 

A(q',t',q,t) = (q'\e- iH ^\q) (6) 



contains all physical information. 

Let us divide the time interval t' — t into n + 1 intervals of equal length 8, 



t'-t 



n+1 



8 is a small quantity in the limit of large n. We can write 



e -iH(t'-t) _ e -iH8 e -iH8 e -iH8 



(7) 



the product of n + 1 factors. 

Again for the sake of notational simplicity we shall assume that 



(8) 



The generic case with linear terms in p ,or ^-dependent coefficients of p 2 , p only present 
technical complications ,but do not alter the results which follow. 

The amplitude eq (6) becomes, by inserting eq (7) and a number of projectors on a 
complete set / dq\q)(q\ = 1 



a = J ^...^^V-^i^^i^i^j)...^!^^) 

On the other hand ,by use of the well known Baker-Haussdorf formula 

e -iHS _ e -i£iS e -iV(q)8 j-j + ff^g 2 )] 



(9) 



(10) 



so that 



and since 



^i +l \e- iH %) = e- lV ^ S (q i+1 \e- i ^\q i ) [l + 0(S 2 )} 



(q i+1 \e- iP ^\ qi ) = J ^(q M \p)<P\9i)e- i ^ = J |V^ +1 -*H^ = 



finally 



(9i + i\*- m \<Ii 



Inserting eq(l 1) into eq (9) 



= lim / 



J SI ££dt 



(ID 



(12) 



In the limit n — > <*>, 8 — > and eq(l 1) becomes exact. Eq(12) is usually written as 



S[q'/,q,t] 



(13) 



The limit in eq(12) defines an integral over an infinite number of variables (a functional 
integral) if it exists. The amplitude eq(6) is an integral over all paths from (q, t) to {q 1 , t') 
and is defined as a limit of a lattice in time ,as the points of the lattice fill densely the 
time interval (t,t'). 

Exponentials in eq(ll) are analytic functions : the construction allows an analytic 
continuation of the amplitude to euclidean time — t — > T = — it If J(t) is an arbitrary 
function in the time interval (/,/') we define 

f' 

S J = S(q,t,t / ) + J dtq(t)J(t) 

Aj(q',t',q,t) = J \®q\ exp [iSj(q' ,t' ,q,t)] (14) 
In the discretized version, 

rt' n 

\ dtq(t)j(t) = J £q(t i )J(t i ) 

Jt i=l 

We define the amplitude 

A(t / ,t l ...t k ,t) = {q / ,t / \T(q(t l )...q(t k ))\q,t} t <t t <t' (15) 
by the same construction which led to eq(13) as a limit of a discretized integral 

A(t',t 1 ...t k ,t) = J[@q]q(t 1 )...q(t k )e- iS ^ 

It follows from the definition of Aj, eq.(14), that 

Mt\ h... t» t) = r sm S " Sj{tk) Aj{q', t\ q, t) (16) 

We shall now compute the amplitude for immaginary times Aj(q' ' ,—iT;q,iT), in the 
limit T — > oo . We shall assume that the external current J is different from zero only 
in the time interval . —T /2 < t < T /2. Then, calling \E n ) a complete set of states with 
definite energy, 

Aj(q',-iT;q,iT)= £ (q' \e~ HT ' 2 \E n ) (E n - Z \E nl ,i^)j(E nl \e- HT / 2 \q) (17) 

F F 

since by assumption from — T to —T/2 and from T /2 to T the evolution is governed by 
H. In the limit of large T 

Aj{q'-iT;q,iT) ~ e^ 7 ¥o W)^(q)(0\ [ [f%]<rW|0> (18) 



|0) is the ground state and E Q its energy. Excited states are exponentially depressed in the 
sum eq( 17). Modulo a multiplicative constant which will be shown to be unimportant,the 
amplitude of eq(16) is then equal to 

Z[J] = (Q\j[®q]e- s W\Q) (19) 

and,as in eq's(15),(16) 

8 n Z[J] 



(0\T(q(t l )...q(t k ))\0) 



8J( h )...8J(t k ) 



(20) 



Z[J] is known as the functional generator of the field correlators (20). 

Since the field correlators fully describe the theory ,knowledge of Z[J] means solution 
of the system. In particular the euclidean Feynman integral uniquely identifies the true 
ground state of the theory. 

Feynman integral quantization is more fundamental than perturbative quantization. 

For QCD the euclidean functional integral will be 

Z[J} = J [S>A pL ][S)w][S)W Sj[A ^ M (21) 

where the role of q is plaied by the gluon fields ,the quark field \]/ and its hermitian 
conjugate \[r. The integral is defined by discretizing a finite volume V of space time to a 
lattice of spacing a and taking first the limit in which the lattice spacing a — > 0, so that 
the lattice densely fills the volume V ,and finally sending V — > °°. 

In QCD both limits exist . By renormalization group arguments the lattice spacing a, 
at a given value of the coupling constant g, is related to the physical momentum scale 
A L as 

a ~ -Lexp(— \) (22) 

where b is minus the inverse of the first coefficient of the beta function of the theory and 
is positive ,because of asymptotic freedom. Eq.(22) will be valid at sufficiently small 
values of g 2 . 

At sufficiently small g 2 a is small in physical units and the granularity of the lattice 
becomes irrelevant. 

Moreover a mass gap exists in the theory (a non zero minimum mass) , which deter- 
mines a finite correlation length ^ . If the lattice is sufficiently large with respect to ^ the 
infinite volume limit is reached. Lattice will be a good approximant to the continuum if 

a <C £, La 

L being the linear size of the lattice. Lattice calculations are in that case an element of 
the sequence eq(12) near to the limit.Lattice will also uniquely identify the ground state. 

The symmetry of the ground state ,which is at the basis of the confinement mechanism 
can safely be studied on the lattice. 



FEYNMAN INTEGRAL AND PERTURBATION THEORY. 



The only known computable functional integral is the gaussian integral,with a la- 
grangean quadratic in the fields, which corresponds to free fields. 
The action can be written as 

5 = X -j dxdy^x^/ix-y)^) (23) 

The resulting equations of motion are 8S/8<p(x) = or / @>^(x — y)<Pj{y) = . The 
matrix ^r. 1 (x — y) is the inverse of the propagator. For the scalar field 

^ ij \ x - y ) = (d 2 + m 2 )8 lJ 8(x-y) 
The discretized version for S Q is 

L ab 

is a symmetric matrix which can be diagonalized by an orthogonal transformation 

T 9- l O = ^ (24) 
The jacobian of the transformation is equal to 1 so that 

Z = J Hdve-^H^h = c(det^) 1 / 2 (25) 
The correlation functions are given by 

(o|r(<p. (xj) . . • |o> = \ J n d( Pi<pH (*i) ■ ■ ■ (26) 

The integral is easy to compute,again by diagonalizing ,giving 

c(det^) 1 /2 £ 

perm 

The sum runs over all possible choices of the pairs i a , L. 

The only non zero correlators are the two point functions. The Hilbert space is Fock 
space. The generic correlator is the product of two point functions. 

In the perturbative approach to quantization the action S is split in a quadratic part 5 
and terms containing higher powers of the fields 

S = S Q + Sj 

The functional integral is then computed by expanding the weight in powers of Sj 

Z= f @(p e - s o- s i= [ $(p e - s ^^—j-S n i ( 27 ) 
J J n n. 



Sj is a polynomial in the fields. The generic correlation function is then given by 



(om<p. ( Xl )... cp in (x n ) \o)=y n^-^H, w • • ■ <p in ( x n) e ^psf (28) 



and again is a gaussian integral which can be computed. 

The result, a sum of products of propagators, is nothing but the well known Wick's 
theorem,and each term corresponds to a Feynman diagram. 

Eq(26) is the interchange of the order of two limits , which is not always legitimate. 

The perturbative expansion can be seen as the functional version of the saddle point 
method for evaluating ordinary integrals. Consider a one variable integral of the form 

rb 

1= / dxexp(if(x))g(x) (29) 

If / and g are analytic in a domain of the complex plane including the real interval {a, b) 
the path of integration can be deformed in that domain to include the points x where the 
phase fix) is stationary,i.e. % = 0. The neighborhood of these points will give the main 
contribution to the integral :far from it the phase factor is rapidly oscillating ,and the 
contribution to the integral negligible. The saddle point method consists in writing / in 
the form 

/= t dde^^f'^gix + d)^^!) (30) 

Ja,C n n - 

where R(8) = f — f(x) — \ f'(x) S 2 (the linear term drops at x ), and the factor multiply- 
ing the exponential is intended as a power series in 8. The integral is then gaussian . If 
more saddle points x i exist / will be the sum of analogous expansions around x i . 

The method works in practice. No systematic control ,however,exists of the approxi- 
mations involved. 

For the Feynman integral the idea is the same, except that now there are infinitely 
many integration variables. The phase is the action S, the integration variables are the 
fields, and a saddle point 

8(p(x) 

is nothing but a classical solution of the equations of motion. 

The point <p = is certainly a saddle point. Around it one can expand the action as 



or i or-! 1 <5 2 S 

S[<p]=S[<p} + 



2 8 <p- Sep j 



StSj + S, (32) 



<p=<p 



For q> = S[<p] = 0, | gf^L 8 t 8j = S and the saddle point formula is identical to the 

perturbative expansion eq (27). 

In general, if different saddle points exist,the expression for the partition function 
becomes 

-5-5. (-1)" 



J>- 5 ^ f[dS]e 2 ^ 5 ^£L_LL 5 « 
if. J „ n. 



and all the classical solutions with finite action contribute. 

Such solutions (if non trivial) are called instantons. In a non abelian gauge theory 
finiteness of the euclidean action means 



d A xTr [G liV (x)G flv (x 



< oo 



or that the field G^ v decreases more rapidly than 1/r 2 as r — > °°. The fields at large 

distance are a pure gauge: = id^U U* with U an element of the gauge group. A 
topological number exists which counts how many times in this mapping the group is 
covered when the point at infinity sweeps the sphere S 3 ■ The topological charge is 



and is an integer. Explicit instanton solutions were found^ which are self dual or 
antiselfdual 

A non zero value of Q implies then a non zero value of G 2 = G^ V G^ V , or that (0|G 2 |0) ^ 
0, a phenomenon which is known as gluon condensation. (G 2 (x)), which must be x 
independent by translation invariance, is known as gluon condensate. A non zero value 
of (G 2 ) is a highly non-perturbative result . Indeed renormalization group dictates for 
(G 2 ) the form of an operator of dimension 4, 

(G 2 ) ~ g-W 2 (33) 

where b > is proportional to the first non zero coefficient of the beta function. An 
expression like eq(33) is non analytic at g 2 = 0, and therefore non expandable in a 
power series of g 2 . A complete classification of the instantonic solutions does not 
exist.Models have been developed based on quasi-solutions,consisting of ensembles of 
single instantons. If their average distance is large compared to their size the ensemble 
is called an instanton gas[0], if the overlapping is important the ensemble is called 
an instanton liquid[||p. Attempts to improve perturbative expansion by adding to the 
perturbative saddle point the saddle points corresponding to instantons have been done 
during the years, starting from the pioneering approach of ref[7]. Although useful to 
describe chiral properties, instantons fail in accounting for confinement of color. 

Instantons were at the basis of the SVZ[|9|] sum rules of QCD: that approach, which 
will be described in next section, is based on the existence of condensates, the origi- 
nal idea being that instantons provide a slowly varying background field for quantum 
fluctuations. 



CONDENSATES. 

The time ordered product of two currents jfi(x),j v (0), e.g. electromagnetic currents,can 
be written at short distances as a sum of local operators times c-number coefficients JlOt]. 

T(j^x)j v (0)) ~ C°. v (x)I + cUx)G a . v G a . v +C*(x)mW+- ■ • (34) 



The operators in eq(34) are ordered by increasing order of dimension in mass. The ex- 
pansion is rigorously valid in perturbation theory, but is assumed to hold also when per- 
turbation theory is not expected to work. Take the vev of eq.(34) and Fourier transform. 
The left side gives 

n MV (<7) = J e^(0\TU^x)j v (0))\0)d 4 x=(q^ v -q 2 gflv )n( q 2 ) (35) 

the tensor structure being dictated by gauge invariance. 
The spectral representation reads 

U(q 2 ) = 11(0) - q 2 f - * ( J° , (36) 



with 



^e+e -^hadrons 

e+e-—fH + iJ.- 



The right side of eq.(34) gives 

(<?M<7v " q 2 g l iv) {C (q 2 ) + C s (q 2 )G 2 + C v (q 2 ) (mxjfyf) + ■ ■ ■} (37) 

The term C Q (q) corresponds to the usual perturbative expansion and is constant modulo 
log's. The other terms, known as "higher twists", are non perturbative in nature. 
The approximate equality 

f°° ds R( 5) 

- q i j T c «2 i iP - {C (q 2 )+C g (q 2 )G 2 + C v (q 2 )(mW) + ■■■} (38) 

JAm^ S S q -tit 

can be exploited by appropriate weighting which emphasises the region where the 
equality is a good approximation on the average:the procedure is known as sum rules[|9p . 

The dispersion integral in eq (38) is evaluated by inserting for R(s) the resonances 
at low s and a constant above some threshold. A nice phenomenology results,relating 
resonance parameters to the condensates, which can be determined consistently giving^] 

G 2 =0.012GeV 4 
(xjfxir) =-0.13GeV 3 

A more recent determination for G 2 isJTT]] 

G 2 = (0.024 ± 0.01 l)GeV 4 



LATTICE DETERMINATION OF THE CONDENSATES. 

Consider the gauge invariant correlator 

^v lP aW = (0|rr{G MV (jc)S c (jc,0)G pff (0)5j(jc,0)} |0) (39) 



where G^ v = £ fl ^-G a ^ v , = £ fl ^A a ^ and X a are the generators of the gauge group in 
the fundamental representation. 



S c (x,0) =P^&xp(iJ* A^(x)dx^ 



is the parallel transport from to x along the path C. Under a gauge transformation U (y) , 
S c (x,0) — > U (x)5 c f/' l '(0), which makes @£ V n a gauge invariant. In what follows we will 

chose for C the straight line connecting and x. ^ v ,pcr can be considered as the split 
point regulator a la Schwinger of the gluon condensate G 2 . The Wilson operator product 
expansion gives indeed 

SfvjHxto C * G 2 (40) 

^ r x^O X 

whence G 2 can be extracted. Poincare' invariance implies [|T2"|] 

%v,pa( x ) = (g f ip8va-8iia8vp)[D(x)+D l (x)]+ (41) 
J f\ x p x p8va x v x p8po x p x o8vp J r x v x o8pp) 
or,choosing the x° axis along x^ 1 one can define [p"5[] 

^|| = ^XXo*=^l+* 2 ^ 

hi 

On the lattice,which is a formulation in terms of parallel transports, the correlators are 
easy to implement: the field strength is the open plaquette, the parallel transport the 
product of elementary links. The correlator is measured as the average on the ensemble 
of the configurations produced by MonteCarlo simulations of the operator: 






c$ L = I 1 > I I L 

-*>VPCT \ N c 

Ppa Ppv 

The figure represents the product of elementary links. The correlators D and D t are 
plotted versus the distance in fm in fig. 1. 
The lines represent a best fit of the form 



D = Ae- x/X + -- I (42) 
\ x \ 

D l = Be ~ X/X + U4 

\Ji 
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FIGURE 1. D (upper curve) and D l vs distance in fm, as determined on lattice[[l3[ [P] ], 



which corresponds to the OPE eq.(40). One can extract A, B and the correlation length 
A . G 2 is proportional, with known coefficients to A + B. The result is 



Quenched[[T3p SU(2) G 2 

Quenched^ ^§ SU(3) G 2 

Full QCD, extrapolated G 2 
to physical q masses [^6]] 



(.33±.01)GeV 4 X 
(.15±.03)GeV 4 X 
(.022±.005)GeV 4 A 



.16 ±.02 fin 
.22±.02fm 
.34±.03fm 



The value of G 2 for full QCD agrees with recent phenomenological determinations [[Tl]] . 
The correlation length is definitely too small in quenched theories, both SU(2) and SU(3) 
to consider the background field as slowly varying, it is slightly larger in full QCD. 
By similar techniques the chiral condensate can be extracted from the correlator[[F7]] 



%{x) =m(\jf(x)S c (x,0)\if(0)) 



which allows the OPE 



C 



%{x) ~—L+C 2 (m\fnir) 



THE STOCHASTIC VACUUM. 



In the Schwinger gauge (x — z Q ) ^A^ (x) = 0, the gauge field A^ can be expressed in terms 
o 



'liv 



A U X ) = J ada(x-z )G£ v (z + a(*-z )) 



(43) 



and any physical quantity can be expressed in terms of correlators of field strengths in 
the vacuum. A cluster expansion is made, and the assumption that higher correlators are 
products of two point correlators, higher clusters being negligible (stochastic vacuum 
model). The two point correlators determined on the lattice are used as an input, and 
several physical quantities can be computed: the heavy qq bound states or high energy 
cross sections. A successfull phenomenology results Jl8|]. 

We have shown that non perturbative effects, like gluon condensate, are important 
in QCD and can be studied on the lattice. In the next lecture we shall investigate what 
symmetry of the vacuum is at the origin of color confinement. 



CONFINEMENT OF COLOR. 

A deconfinement transition is observed in lattice simulations of QCD. 

The static thermodynamics of a field system at temperature T is described by the 
partition function 

Z = tr[e- H ' T ] 

It is a known theorem that Z is given by a Feynman euclidean integral, extending on the 
3-dimensional physical space and on the time interval (0, l/T), with periodic boundary 
conditions for bosons,antiperiodic for fermions. 

On the lattice Z is obtained by simulating the theory on a lattice N$ x N T ,with 
N s ^> N T . In terms of the lattice spacing a the temperature is then given by T = 1 /N T a. 
The only free parameter in the simulation is the gauge coupling constant g 2 or better the 
parameter /3 = 2N c /g 2 . In terms of /3 the lattice spacing can be written as [see eq(22)] 

a = x exp(--) 

or 

T = — = — exp(-) (44) 
N T a N T b 

Low temperature corresponds to small /3 (strong coupling), high temperature to large /3 
(weak coupling). 

This is the opposite to what happens in ordinary spin systems, where the temperature 
plays the role of coupling constant, and is due to asymptotic freedom. Confinement, 
for pure gauge theories, is related to the Polyakov line, which is the parallel transport 
along the time (temperature) axis from to N, closed by periodic boundary conditions. 
The vev of the Polyakov line (L) can be interpreted as exp(— fx q ), \x q being the chemical 
potential of an isolated quark. In the confined phase ji q = °o, (L) = 0, while one can have 
(L) 7^ in the deconfined phase. (L) can be called an order parameter for confinement, 
the symmetry being Z N . What is usually determined on the lattice [|19|] is the correlator 
of two Polyakov lines, which by cluster property behaves at large distances as 

G(d) = (L(x)L(y)) ~Aexp(-odaN T ) + \(L)\ 2 (45) 

with d = \x — y\. 



The energy of two static quarks at distance d is related to G(d) by the equation 



V(d) = -^-log(G(d)) (46) 
aN T 

What is found on the lattice is that a temperature T c exists such that 

T <T C |(L)| = V(d)- d ^ oo od (Confinement) 
T > T c | (L) j / V(d) const. (Deconfinement) 

For quenched SU (2) one finds T c / yfo ~ 0.7. 

A finite size scaling analysis of G at different spatial sizes N s , allows to extract the 
critical index v The result is v ~ 0.62. The phase transition belongs to the same class of 
universality as the 3d ising model,as expected fl20|]. For quenched SU (3) the transition is 
found to be weak first order, v = 1/3 and T c / \fd = 0.65 which means T c = 270 Mev if 
the conventional value y/o = 425 MeV is assumed for the string tension. 

In the presence of quarks (L) cannot be an order parameter, since the symmetry Z N 
is explicitely broken by the very presence of the quarks. At m q = 0, the chiral limit, an 
order parameter is the chiral condensate (fyy): however in reality the chiral symmetry 
is explicitely broken by quark masses. This variety of order parameters contrasts with 
the ideas of the limit N c — > °°. 

The idea is that the number of colors N c can be considered as a parameter of the 
theory. In the limit N c — > °° with g N c = A fixed, a theory is defined , which differs little 
from the theory at finite N c , say N c = 3. 

In this spirit the mechanism of confinement should be fixed by the limiting theory 
,and be essentially the same at all values of /V c . 

Quark loops being non leading in the expansion ,the mechanism of confinement 
should be the same in full QCD and in quenched theory. 

The symmetry of the vacuum which is behind confinement, as explained in sect (1) 
should be /V c - independent. In this respect the existing situation for the order parameters 
looks rather confusing. 

Fig. 2 shows schematically the phase dyagram for 2-flavor QCD as a function of the 
quark masses. The border between confined and deconfined phases is defined as the 
simultaneous maximum of various susceptibilities[|TJ. >From m q = <=o to m q = 3 GeV 
the transition is first order, and the Polyakov line works as an order parameter as in 
quenched QCD. At m q = there are arguments that the transition is second order[^2|]. 
In the central region none of the susceptibilities which have been determined (ty/ty/, 
susceptibility, derivative of the specific heat) goes large at the border between the two 
regions as the volume goes to infinity, and the conclusion is that maybe there is no 
transition, but only a crossover. 

Still, if the deconfining transition is a change of symmetry,as argued in sect.l, there 
must exist an order(or disorder) parameter which distinguishes confined from decon- 
fined, possibly the same as in quenched theory,in the line of /V c — > °° ideas. 



300 




m =3GeV m =00 

q 4 



FIGURE 2. A schematic phase dyagram for 2 flavour QCD. Shaded area is confined, upper area 
deconfined. 

DUALITY 

The symmetry responsible for confinement has to be a symmetry of the (strong coupling) 
confining phase, which, in the langauge of statistical mechanics, is disordered. The 
natural question is: what can be the symmetry of a disordered phase? The answer is 
in a known concept in statistical mechanics, namely duality [[23]]. d dimensional systems 
admitting non local configurations with non trivial topology in d — 1 dimensions, can 
have two complementary descriptions. 

One description in terms of the ordinary local fields 0,in which topological excitations 
/i are non local, and (0) are the order parameters; this description is convenient in the 
weak coupling regime g <C 1 (ordered phase). 

The other one ,dual description,in which the excitations \x become local fields ,their 
vev's (n) are the order parameters: the coupling constant g D in this description is related 
to the original g as g D ~ l/g. Duality maps the strong coupling regime of the direct 



m q = 



description into the weak coupling regime of the dual. 

A number of systems admitting dual descriptions are known in the literature.The 
prototype system is the 2d Ising model [^4]], in which the dual excitations are kinks, and 
the dual description is again a 2d ising model with g D ~ 1/g; the N=2 SUSY QCD of 



refQ25|], in which dual excitations are monopoles; the 3d XY model Q26|] where the dual 
excitations are vortices, and the Heisenberg magnet in 3d, where the dual excitations are 
2d Weiss domains [|27|]. In QCD the dual excitations are not exactly known, but important 
information on their properties exists as we shall see below. 
Two main candidates were originally proposed: 



a) monopoles pq, |29| ]. Their condensation in the confining vacuum generates mag- 
netic superconductivity,in the same way as the condensation of Cooper pairs gen- 
erates ordinary superconductivity. 

The chromoelectric field between a qq pair is channeled into electric Abrikosov 
flux tubes, so that energy is proportional to the distance. 

b) Z N vortices [|3TJ 



In what follows we shall analyze the option a). As for vortices we refer to[[3T| |32p and to 
references therein. 



BASIC SUPERCONDUCTIVITY^. 

Before looking for dual (magnetic) superconductivity, we recall the basic concepts of 
ordinary superconductivity. The Landau-Ginzburg[|34j] density of free energy (effective 
lagrangean) is determined by arguments of symmetry and of scale. 

1 1 tn 2 X 

% = -4 W + ^{D^YD^ + ^-0*0 - -(0» 2 + irrelevant terms (47) 

is a charged scalar field describing Cooper pairs , = <9„ — iqA^ is the covariant 
derivative and the irrelevant terms are those of dimension higher than 4. m 2 = m 2 {T) 
and X = X{T) depend on the temperature. 

For T < T c m 2 > and the potential has a mexican hat shape as a function of R((j)) an 
of 7m(0). 

For T > T c m 2 < and the minimum of the potential is at = 0. If we parametrize 
^ = pe' qd , then = — d^6 is a gauge invariant quantity: indeed under a gauge 
transformation A^ — ► A^ + d^A, — > 6 + A. In terms of the new variables 

^ = ~\f, v F, v + ^fA,A, + 1 -(d,p) 2 + ^p 2 - \p" (48) 

A static homogeneous solution is p = p, a constant. The effective lagrangean for the 
photon is in that case 

^ = - l -F ilv F ilv +q 2 p 2 A il A il 

The equation of motion is 

V AH + q 2 p 2 A = (49) 



Taking the curl of both sides of eq(49) 



V 2 H - q 2 p 2 H = (50) 

The second term in eq (49) is a gauge invariant current,and is known as London current. 

A non zero current with zero electric field means zero resistivity or superconductivity. 

Eq(50) means that the magnetic field has a finite penetration depth in the supercon- 
ductor , l/qp, and this is nothing but Meissner effect. If m^> qp the superconductor 
is named Type EL In that case, when trying to introduce a magnetic field in the bulk, a 
penetration through separated Abrikosov flux tubes is energetically favoured. 

Outside the tube A = and § c Adx = 0, C being a path which encercles the flux tube. 

Since A = A — V<j>, this means 

/ Adx= / Vdx= (51) 

JC JC q 

The magnetic flux in the tube is quantized (Dirac quantization condition). A flux tube is 
to all effects a monopole antimonopole pair sitting at the ends with energy proportional 
to the length. 

The order parameter of the system is <I> = qp, the vev of a charged field. 



MONOPOLES IN QCD. 

Monopoles are always abelian[R5]]: the magnetic monopole term in the multipole expan- 



sion of the field produced at large distances by any hadronic matter distribution ,obeys 
abelian field equations, and can always be reduced by a gauge transformation to the form 
(0, 0, a/r) in polar coordinates, with 

a 9 = 2(1 -cose) (52) 

and the Dirac string along the north pole. The N c x N c matrix Q obeys the condition 
e ig4xQ _ J ^ implying that,in the representation in which Q is diagonal, it has integer or 
half integer eigenvalues in units l/g. A monopole is identified by a constant diagonal 
matrix of the algebra with integer or half integer matrix elements: there exist iV c — 1 
independent monopole charges. 

The same result follows from the procedure known as abelian projection[p8|]. We 
shall illustrate it for the case of SU(2) gauge group: generalization to arbitrary N c is 
straightforward. Let be any field in the adjoint representation. Define = 0(x)/|0(jc)|, 
a direction in color space, and 

Fp V = f <5 MV - -$ (D^AD V $) (53) 

§ 

Both terms in the expression (53) are color singlets and gauge invariant.The combination 
is chosen in such a way that bilinear terms A^a v and A^dyty cancel. Indeed, by explicit 



computation 

F liV = $(d ll A v -d v A ll )--$(d ll $Ad v $) (54) 

A gauge transformation which brings $ to a constant, e.g. to (0, 0, 1) is called an abelian 
projection on $(x). In that gauge indeed 

F^ v = d^A v )-d v (^) (55) 

is an abelian field. If, in the usual notation F* v = £^ V p a F p a, a magnetic current can be 
defined as 

This current is identically conserved, because of the antisymmetry of F* v 

d v j v =0 (57) 

and identifies a U(l) magnetic symmetry of the system. 

In the usual continuum formulation j v is zero (Bianchi identities). 
In a compact formulation, like lattice, j v can be non zero. 

The U(l) magnetic symmetry can be either realized a la Wigner, and then the Hilbert 
space consists of superselected sectors with definite magnetic charge, or can be Higgs 
broken. In the first case the vev (fl) = (0|jix|0) of any operator \x carrying magnetic 
charge is strictly zero. In the second case there exists some \x such that (/i) ^ 0. 

The effective action for (/i) will then be of the form eq(47) and the system behaves 
as a dual superconductor. The operator /i is a disorder operator 

(/i) 7^ in the disordered,confining phase 
(jii) = in the ordered, deconfined phase 

A magnetically charged operator ji can be constructed,which is magnetically charged 



in any given abelian projection[j36|, p7[ ]. The basic principle to construct \i is the well 
known formula for translations 

e ipa \x) = \x + a) (58) 

In field theory the field 0(y) is the analog of x. Its conjugate momentum n(y) the analog 
of p, and the operator 



ij (Pyliiy^ix-y) 



(59) 



H(xj) = exp 

gives 

ii(x,t)\ii>(y,t)) = \^(s,t) (60) 

i.e. it adds the classical configuration (j)(x — y) to any field configuration. 

For monopoles the field is A(x,t), the conjugate momentum <9qA(x,?) = E(x,t), and 
the classical configuration to be added is the vector potential generated at the point y by 
a monopole sitting at the point x , e.g. 



where the Dirac string has been put along the n axis. The definition (60) has to be adapted 
to compact formulation .The details are published in ref P7|] . The result is of the form 

where Z(j8) is the partition function 

Z(J8) = |[«A][^[^]«ap(-J8S £ ) 

and 

Z(/3) = f[®A][9ir][®¥]exp{-ps E ) 



S E differs from S E by the introduction of a dislocation on the slice x° = t. AS is 
proportional to the spatial volume, and therefore fluctuates as N^ 2 , so that Z fluctuates 
as exp(A^/ 2 ). It proves convenient to define instead 

p(j8) = ^log(M) = (5) 5 -(5 + A5) 5+A5 (62) 
which is easier to measure.In terms of p 

(ju)=exp / p(x)dT (63) 
Fig. 3 shows schematically the shape of (ju) as a function of /3, and that of p . 




FIGURE 3. Typical dependence on j3 of ji and p. 



For /3 < /3 C (confinement) (ju) ^ 0, for /3 > j8 c 

(M)~(/3-/3 c ) 5 = (l 



</x> = 0, at j8 - j8 c 

7/ 



For finite lattices (fl) cannot be strictly zero above /3 C ,because it is analytic in /3 and 
therefore if it were zero on a line it would be identically zero. Only in the infinite 
volume limit singularities can develop^] and can be identically zero above /3 C . The 
limit N s — > °° can be studied by finite size scaling techniques. 

I In the region of low /3, at N s larger than any physical length scale p becomes N s 
independent and finite, and from eq(63) one concludes that (/x) ^ (fig. 4). 
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FIGURE 4. Limit at large N s of p for j3 < J3 C , 51/(2) gauge group. 

II .In the region /3 > /3 C - the dependence of p on N s can be studied . It is found (see 
fig.5) 

p = -kN s + k' k>0 (64) 

In the limit N s — > °°, (jii) is strictly zero. A direct determination of (/i) would give 
zero within large errors. 

Ill In the vicinity of j8 c by dimensional arguments 

(V) = ? S <S>(^J) (65) 

where ^ is the correlation length in units of lattice spacings. ^ diverges at the critical 
point with an index v 



T 



~ %~ v x = 1 - ^ (66) 
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FIGURE 5. p vs N s for j3 > j3 c . 



If 4 » 1 then l/§ ~0and 

(M)~T 5 4»(0,^)^ ~ T- v (67) 
One can trade 7Y 5 / ^ with v x and 

p = -- +n¥ v &(0,nV v t) (68) 

or 

p/Ny* = -* + <6r(0,N}/Vx) (69) 

A scaling law results: p /Nl' v is a universal function of N\I V T, independent of N s . Fig. 6 
shows a typical dependence of p on A/j. Fig.7 shows the quality of scaling for SU(2). A 
best square fit to the data allows to determine /3 C , v, 8. 
For pure gauge theory one obtains p5|] 

SU(2) V = .62(l) 5 = 0.20(3) 
SU(3) v = .33(l) 8 = .50(3) 

The result is independent on the choice of the field used to perform the abelian 
projection. Confining vacuum is a dual superconductor in all abelian projections[]3~8|]. 
In all projections magnetic symmetry is Wigner in the deconfined phase. 
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FIGURE 6. Volume dependence of p, pure gauge SU (2), deconfinement transition. 



Whatever the dual excitations of QCD are, they must be magnetically charged in all 
the abelian projections. 

In QCD with dynamical quarks the same operator (fi) can be defined ,and the corre- 
sponding p = 4j log((/i}) can be studied Results already exist that (jti) ^ in the 

confining phase, and (ju) = in the infinite volume limit. . The finite size scaling in the 
vicinity of the transition is on the way. 

Now another scale, the quark mass enters, which was absent in the quenched case. 
Eq(65) -(67) become 

(M)^T 5 4>(i J-,m^) ~ T 8 <$>{0,N-s l l v T,m q N^) (70) 
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FIGURE 7. Scaling, eq.(69). ST/ (2) pure gauge. 



The index y is known. In order to determine V one can choose values of N s such that 
m q^l' /r =const., and test the scaling law 

p/N l J v ~<$>(N 1 J v t) (71) 

to determine v and to investigate the nature of the transition. This will complete the 
analysis. The existing data already show,however,that dual superconductivity is the 
mechanism of confinement also in the presence of dynamical quarks, in agreement with 
the ideas of N c — > °°. Fig. 8 shows the negative peak of p at the same temperature where 
the chiral condensate drops to zeropO[|. 



CONCLUDING REMARKS. 

Large distance QCD is an intrinsecally non perturbative system.Lattice provides a tool 
to investigate this regime. 

Much progress has been done towards the understanding of confinement. The mech- 
anism of confinement is definitely dual superconductivity both for pure gauge and full 
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FIGURE 8. 2 flavour QCD. Open circles represent the chiral condensate, with the scale axis on the left. 
Full circles represent p, with the scale on the right. The peak of p coincides with the drop of (xj/ij/). 



QCD. The dual excitations which condense to produce confinement are not yet identi- 
fied: what is known is that they are magnetically charged in all abelian projections. 
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